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We investigate the role of lattice polarization in determination of induced carrier density at the n-type in- 
terface of LaA103 overlayer on SrTiO;i (001) by carrying out density-functional-theory calculations. When no 
oxygen vacancy or defect is present, the magnitude of polarization screening in the LaAlC>3 layers is found to 
be correlated with the carrier charge induced at the interface. For the interfaces with a few LaAlC>3 layers, the 
induced charge carrier is compensated by the electrostatic screening and consequently its density remains far 
less than 0.5 electrons per unit cell. 

PACS numbers: 73.20.-r, 79.60.Jv, 77.22.Ej, 73.21.-b 



The observation of a high mobility electron gas in the n- 
type (LaO)/(Ti02) interface between two band-gap insulators 
LaAlOa (LAO) and SrTiOa (STO) has generated intense re- 
search activities toward its potential device applications as 
well as its physical mechanism 1 1]. The electronic reconstruc- 
tion at the interface has been suggested as a way to avoid the 
polarization catastrophe yl] which may arise from the alter- 
nating stack of positively charged LaO and negatively charged 
AIO2 layers on top of the TiC>2 termination of STO substrate. 

While the reports of ferromagnetism y[ and superconduc- 
tivity J3l in the LAO/STO interface have boosted up research 
interest in the mechanism of conductivity and the dimension- 
ality of the induced charge carrier at the interface, there still 
remain controversies on the origin and nature of the interface 
electron gas Several experiments have demonstrated 

that oxygen vacancies in the STO layer are responsible for the 
high carrier density, which depends more on the film growth 
and annealing conditions J3, d, 1st] . Apart from the oxygen- 
vacancy-generated carriers in samples grown in oxygen-poor 
conditions, the carrier density of the n-type interfaces with 
perfect stoichiometry poses another puzzle that a common 
lower limit in the carrier density at ~ 10 13 , i.e., 0.03 electrons 
per unit-cell (u.c), is order-of-magnitude less than that of 0.5 
electrons per u.c. expected from the electronic reconstruction 
mechanism ifioll . 

On the other hand, recent density-functional-theory (DFT) 
calculations pointed out the importance of polar distortions as 
a source of dielectric screening in the LAO/STO heterostruc- 
ture ll ll Il2ll . Without the lattice relaxation, even a single 
layer of LAO on STO(001) with an ideal structure would be- 
come metallic. Further the insulator-to-metal transition can be 
driven by an external electrical field [l3ll . The thickness 
dependence of carrier density also suggested a possible role 
of electrostatic screening in the LAO/STO interface 11411 . 

Here we show that the detailed balance between the lattice 
polarization and the charge transfer is important in determi- 
nation of the carrier density at the n-type interface of LAO 
overlayers on STO(001) with perfect stoichiometry. Our DFT 
calculations demonstrate that the lattice polarization of the 
LAO overlayer is correlated with the carrier charge induced 



at the interface in terms of the overlayer thickness. When the 
LAO overlayers are over a critical thickness, the charge trans- 
fer from the LAO surface to the interface is compensated by 
the electrostatic screening due to the polarization distortions 
across the LAO layers. 

We carried out DFT calculations by using the Vienna ab 
initio Simulation Packa ge ( VASP) IU5I1 within a generalized 
gradient approximation 11611 together with the projector aug- 
mented wave pseudopotentials Jl5l 17 1 and the cut-off en- 
ergy of 400 eV for the plane wave basis. We modeled the 
LAO/STO interface by a slab consisting of one-to-ten LAO 
layers on top of seven STO(001) layers and a vacuum region 
of 14 A along the c-axis in a supercell geometry. Dipole cor- 
rections were used to correct the errors of electrostatic poten- 
tial, forces, and total energy caused by periodic boundary con- 
dition lfl9ll . The in-plane unit cell for calculations was taken 
as (lxl) because the rotation distortion of TiOg octahedron 
is known to be negligible for small carrier doping 1 12. 20ll . 
The in-plane lattice constant of the slab was constrained at 
the calculated equilibrium lattice constant a = 3.942 A of 
the STO substrate. The c-axis coordinates of atomic positions 
were fully relaxed with forces less than O.OleV/A except for 
the atoms in the bottom three layers of STO, which were fixed 
in their bulk structure. 

In order to determine accurately the amount and distri- 
bution of carrier density as a function of the layer thick- 
ness, we had to devise a proper k-point sampling for the 
LAO/STO slab calculations. For insulating states, the k-point 
grid of(12xl2xl) was found to be good enough for the self- 
consistency iterations. For the case with metallic charge carri- 
ers, however, it was required to introduce an extra dense grid 
set near T and M points because the accurate evaluation of 
the tiny pocket sizes of Fermi surfaces is critical for the cor- 
rect description of the charge carrier as well as the electrostatic 
potential across the interface. 

One of the most prominent features in the calculated re- 
laxation geometry of the LAO/STO slab is the presence of a 
large polar distortion in the LaO and AIO2 layers but a neg- 
ligible distortion in the STO side B12I1 . Because of the huge 
dipole field produced by the alternating stack of (LaO) + and 



2 



(AIO2) - charged layers, all the LaO and AIO2 layers become 
buckled except the surface AIO2 layer which exhibits only a 
small uniform relaxation. The relative displacement between 
La and O ions within the LaO layers varies from 0.1 A to 0.35 
A, depending on the LAO thickness as well as the location of 
LaO layers, while the relative displacements between Al and 
O ions within the AIO2 layers remain less than 0.15 A. 

The screening by the polar distortion is critical in under- 
standing the interface electronic structure. We obtained an 
interesting relaxation dependence of the LAO/STO interface 
character in the overlayer structure, which is similar to the 
results by Pentcheva and Pickett ll 211 but apparently quite dif- 
ferent from that of (LAO) m /(STO)„ superlattices JUl. Re- 
gardless of their thickness, all the unrelaxed overlayer struc- 
tures were found to have a metallic interface. After the relax- 
ation, however, the LAO/STO interface remains insulating up 
to four layers of LAO. In the systems with one-to-four LAO 
layers, the lattice polar distortion is favored energetically over 
the charge transfer between the surface and the interface. 
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FIG. 1: (Color online) Layer projected density-of-states (pDOS) of 
(a) (LAO) 3 /(STO) 7 and (b) (LAO) 6 /(STO) 7 . Dashed lines are drawn 
by connecting band edges to represent the layer-dependent potential 
profile. 

Figure Q] shows the layer-projected density-of-states 
(pDOS) of two representative systems with fully relaxed LAO 
layers: (a) (LAO) 3 /(STO) 7 and (b) (LAO) 6 /(STO) 7 . The 
pDOS of the relaxed (LAO)3/(STO) 7 exhibits an insulating 
band structure where the valence band maximum (VBM) of 
the LAO layers is set at the LAO surface layer. The dipole 
field in the LAO layers shifts the LAO band edges toward the 
higher energy relative to those of STO. The potential gradient 
is about 0.6 eV per layer for the relaxed structure, whereas it 
would be about 1.3 eV for the unrelaxed structure without po- 
lar distortions. It is noted that the VBM of the LAO surface 
layers, i.e., LAO-3 of Fig. [Da) and LAO-6 of Fig. Qlb) de- 
viates significantly from the extrapolated (dashed) lines con- 



necting the VBM's of its sub-surface layers. The lowering 
of the VBM of the LAO surface layer is attributed to its lo- 
cal environment being close to the bulk AI2O3 with a band 
gap larger than that of LAO. In addition to the surface ef- 
fect, there also exists a band offset at the interface. There is a 
small but clear offset of about 0.5 eV between the LAO-1 and 
STO-1 layers, while the conduction band minimum (CBM) 
of the LAO-1 interface layer smears due to the hybridization 
between LaO and Ti02 layers. 

The band edge shift in the LAO side reflects the change of 
electrostatic potential across the LAO overlayers. In Fig. Q] 
the layer-dependent potential profiles are marked by dashed 
lines. The gradient of the edge shift in (a) (LAO) 3 /(STO) 7 
is higher than that in (b) (LAO)e/(STO) 7 . From our calcu- 
lations, in the system with m = 4 in (LAO) m /(STO) 7 , the 
VBM of the surface LAO layer barely touches the CBM of 
the interface STO layer. When we have more than four LAO 
layers, the LAO/STO interface becomes metallic. As shown 
in Fig.QIb) of the pDOS of the relaxed (LAO) 6 /(STO) 7 struc- 
ture, the charge carrier at the LAO/STO interface is induced as 
a result of charge transfer from the LAO surface. Despite that 
the potential gradient is reduced significantly by the polariza- 
tion screening, the total potential shift across the LAO layers 
exceeds the STO band gap when the overlayer thickness goes 
over a critical value of 4 LAO layers. If the VBM of the sur- 
face LAO layer lies above the CBM of the interface STO layer, 
a charge transfer can occur between the LAO surface and the 
LAO/STO interface. From the results, we found that the total 
shift of the VBM is always delimited by the STO band gap for 
all m > 4. It implies that the potential shift is somehow com- 
pensated by the charge transfer. The polar distortion occurs 
in response to a local electric field, i.e., the potential gradient 
in the LAO layers, which in turn depends on the amount of 
charge transfer. Therefore, in order to determine the amount 
of charge transfer, it is necessary to understand the relation 
between the lattice polarization and the charge transfer. 

In order to quantify the thickness dependence of the lat- 
tice polarization, we calculated the layer-by-layer polarization 
over the LAO layers as a functional of the LAO layer thick- 
ness. Polarization Pg of the ^-th layer can be estimated by 
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where Z* H and Sua are the Born effective charge and the dis- 
placement of the i-th ion in the £-th layer respectively, and ilg 
is the volume of the LAO unit cell in the ^-th layer. Here we 
took the values of the Born effective charges from Ref. |22. 
Since Eq. (fTJ can be used for the either La or Al cation- 
centered cells, we averaged out the contributions of the oxy- 
gens in the planes bounding the elementary unit cells 12311 . 

Figure [2] illustrates the profile of the layer-by-layer polar- 
ization P( over the LAO layers for (LAO) m /(STO) 7 (m = 
1 ~ 7). In all cases, it is noted that the polarizations of the in- 
terface and surface layers are significantly different from their 
average values. Although a reduction of the interface polar- 
ization may arise from both structural relaxation and charge 
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FIG. 2: (Color online) Calculated layer-by-layer polarization Pg over 
the LAO layers for (LAO) m /(STO) 7 (m = 1 ~ 7). Labels "La" and 
"Al" correspond to the LaO and AIO2 layers respectively. 



screening, its origin may require a further investigation but is 
expected to be different from that of the ferroelectric insulator 
superlattices 124ll . Aside from the polarizations at the inter- 
face and surface layers, the average values of Pg for n = 2, 3, 
and 4, denoted by a dashed line in Fig. [2] are roughly constant 
at about 0.03 e/A 2 . From m = 4 to 7, The average value of Pg 
drops quickly as the layer thickness increases. This decrease 
of Pg is related to the reduction of the potential gradient in the 
LAO layers, which is delimited by the STO band gap and the 
charge transfer between the surface and the interface. 

For a given thickness t for m > 4, since the potential shift 
in the LAO layers is delimited by the STO band gap, the local 



electric field Ed can be approximated by Ed = V^ STO / (t—to) 
where V~ is a potential change corresponding to the STO 
band gap, and to is an effective parameter corresponding to 
the dead layer contribution 12511 . As the layer polarization Pg 
is proportional to the local field Ed, the inverse of the average 
polarization Pg should satisfy the relation: 
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where \ is an average polarizability of the LAO layers. The 
t-dependence of Pg in Eq. (O is clearly demonstrated in the 
plot of 1/ (Pg) vs. t as shown in Fig.|3ja). For m > 4, 1/ '(Pg) 
increases linearly in t, while it remains constant below m = 4 
except for the m = 1 case. The deviation of the m = 1 polar- 
ization is not surprising because the single LAO layer consists 
of both surface and interface layers. From the least square fit 
to the calculated data, we could determine the effective pa- 
rameter to by to = 0.053co where cq is the unit-cell thickness 
of the LAO layer. 

For the metallic interface with m > 4, the potential shift, 
generated by the dipole field arising from the alternating stack 
of (LaO) + and (A102)~ charged layers, is screened by both 
lattice polarization and charge transfer. The relation between 
the interface (or surface) charge a c and the local electric field 
Ed can be described in a simple electrostatic picture by the 
electric displacement field across the LAO overlayer: er c = 
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FIG. 3: (Color online) The LAO thickness dependence of (a) the 
average inverse polarization 1/ (Pg) and (b) the charge carrier density 

n c . 



eEd = €oEd + P- If there is no free charge carrier at the 
interface, an average displacement field should correspond to 
<j c = 1 / 2 as expected from the polarization catastrophe theory 
Jll]. When a charge transfer of n c occurs, however, the charge 
screening should give rise to a c = 1/2— n c . By combining the 
t-dependence of Ed and Pg in Eq. (|2), we obtain an expression 
for the t-dependence of the induced charge carrier n c at the 
interface for m > 4: 
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where A 

The calculated charge carriers density n c at the LAO/STO 
are shown in Fig. Ob). For m < 4, the insulating interface 
has no charge carrier. At m = 4, where the VBM of LAO 
touches the CBM of STO, the carrier density n c is ~ 10~ 3 . 
For m > 4, n c increases monotonically as the thickness of the 
LAO overlayer grows. For the thin layers with m < 12, the 
carrier densities remain far less than 0.5 e/u.c. From Eq. (0, 
it is obvious that n c = 0.5 e/u.c. can be achieved only when 
the LAO layers become extremely thick. By anchoring the 
coefficient A to the value of n c at m = 4, we obtained A = 
1.97. The fitting result based on the electrostatic model of 
Eq. <£3j is displayed as a dashed line in Fig. FJJb). Despite 
of some discrepancy for the thin layers close to m = 4, the 
electrostatic model is found to be in a reasonable agreement 
with the DFT results. Considering that the electrostatic model 
is simplistic and may not be valid in the scale of a few atomic 
layers, this agreement is rather remarkable. 

Although our electrostatic model explains the overall be- 
havior of n c without oxygen vacancy, the induced charge car- 
rier still depends on the detailed electronic structure near the 
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FIG. 4: (Color online) Layer-resolved carrier densities in STO for 
the metallic interfaces with the varying thickness of LAO oveiiayers. 



interface. The charge carrier distribution at the interface is 
closely related to the band bending, i.e., the CBM edge shift 
in the STO substrate. When there is no charge carrier, i.e., 
n c = 0, for the insulating interface with m < 4, no band 
bending occurs in the STO side as shown in Fig.[Tfa). For the 
metallic interface with m > 4, on the other hand, the degree of 
the band bending depends on the amount of the charge carrier. 
In Fig. [4] we present the layer-resolved carrier distribution as 
a function of the LAO layer thickness. The induced charge 
carriers are found to be confined within 5 layers close to the 
LAO/STO interface, while the distribution for m = 5 and 6 
is even more localized within less than 3 layers. This result 
is in good agreement with a recent X-ray photo-electron spec- 
troscopy (XPS) experiment, which reported that the carriers 
far less than 0.5 e/u.c. were found to be confined in one or a 
few layers of STO for thin LAO film without oxygen vacancy 
& 

In conclusion, we explored a microscopic picture for the 
screening mechanism at the interface of the LAO overlayer on 
STO(001) by carrying out the DFT calculations, and demon- 
strated that the lattice polarization of the LAO layers is corre- 
lated with the carrier charge induced at the interface in terms 
of the overlayer thickness. Here we suggest an electrostatic 
model for the description of the charge carrier induced at the 
interface, where the detailed balance between the lattice polar- 
ization and the charge transfer plays a primary role when no 
oxygen vacancy or defect is present in the system. Although 
our discussion is restricted to the perfect lattice without de- 
fect, the electrostatic screening mechanism across the LAO 
overlayer should be considered as an alternative to the elec- 
tronic reconstruction mechanism. We hope that our findings 
contribute to resolving the controversies on the origin and na- 
ture of the interface electron gas. In this picture, it is natural 
to observe a common lower limit of 0.03 electrons per u.c. 
for the systems with 4-to-5 LAO layers. In addition, one can 
consider the induced carrier density controlled by an external 
field, which affect the lattice polarization over the LAO layers. 
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